Transport spectroscopy of a single dopant in a gated silicon nanowire 
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We report on spectroscopy of a single dopant atom in silicon by resonant tunneling between 
source and drain of a gated nanowire etched from silicon on insulator. The electronic states of this 
dopant isolated in the channel appear as resonances in the low temperature conductance at energies 
below the conduction band edge. We observe the two possible charge states successively occupied 
by spin-up and spin-down electrons under magnetic field. The first resonance is consistent with the 
binding energy of the neutral D° state of an arsenic donor. The second resonance shows a reduced 
charging energy due to the electrostatic coupling of the charged D~ state with electrodes. Excited 
states and Zeeman splitting under magnetic field present large energies potentially useful to build 
atomic scale devices. 

PACS numbers: 73.21, 73.23, 73.63 



Dopant atoms are essential in semiconductor technol- 
ogy since they provide extrinsic charges necessary to cre- 
ate devices such as diodes and transistors. Nowadays the 
size of these electronic devices can be made so small than 
the discreteness of doping can influence their electrical 
properties [l[ . On the other hand it may be an important 
breakthrough if a dopant could be used as the functional 
part of a device instead of just providing charges. As 
an example, dopant-based spin qubits in silicon are pos- 
sible candidates for quantum computation 0, [3| thanks 
to their longer spin coherence time [3| as compared to 
two-dimensional quantum dots defined by top gates in 
III/V heterostructures Although dopants are well 

known in bulk semiconductors, specific questions arise 
in the context of nano-scale devices like the reduced life- 
time of the two-electron state under electric field involved 
by read-out schemes of spin qubits [7| . The aim of this 
work is to study the electronic states of single dopants in 
gated silicon nanostructures to bring information useful 
for these issues. 

Electron tunneling through isolated impurities has 
been observed previously in two-terminal devices such 
as GaAs/AlGaAs double barrier heterostructures @, Q. 
Here we present experimental results on electron trans- 
port through the localized states of individual n-type 
dopants in silicon nanowires. In contrast to previous 
studies, our devices have a three-terminal configuration 
with source, drain, and gate electrodes allowing a de- 
tailed investigation of charge, orbital, and spin states. 
In particular, we observe both the neutral D° and nega- 
tively charged D~ states, and compare their binding en- 
ergy with the case of bulk dopants. This work provides a 
quantitative description of the electronic properties of a 
single dopant connected to electrodes in a gated nanos- 
tructure. This is the first transport experiment measur- 
ing the charge states of a real atomic system with a 1/r 
attracting Coulomb potential, thus very different from 



quantum dots with harmonic potentials. 

Our devices are 60 nm tall crystalline silicon wires 
(fins) with large contacts patterned by 193 nm optical 
lithography and dry etching from Silicon-On-Insulator. 
After a boron channel implantation, a 100 nm poly- 
crystalline silicon was deposited on top of a nitrided oxide 
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FIG. 1: (color online), (a) SEM images of typical devices, the 
scale bar is 200 nm. (b) Cross-sections showing n- and p-type 
regions. The channel forms at the gate oxide with enhance- 
ment at the top corners, (c) Height of the channel barrier 
measured at high temperatures (squares). The band edge Ec 
reaches the Fermi level Ef at 400 mV. Inset: Low tempera- 
ture differential conductance showing resonances below this 
threshold (al to b2) and above (cl and c2). 
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(1.4 nm equivalent oxide thickness), then received a phos- 
phorus (P) implant as pre-doping, and was patterned us- 
ing an oxide hard-mask to form a narrow gate (Fig.[T]a,b). 
Next, we used high- angle arsenic (As) implantations as 
source/drain extensions, while the channel was protected 
by the gate and 50 nm wide nitride spacers and remains 
p-type. Finally, As and P implants and a NiSi metallic 
silicide are used to complete the source/drain electrodes. 
Conventional operation of this n-p-n field effect transistor 
is to apply a positive gate voltage to create an inversion 
in the channel and allow a current to flow. At low tem- 
perature, the sub-threshold conductance shows series of 
conductance peaks due to electron tunneling through a 
potential well in the conduction band [lOj . This well lo- 
cated in the channel is separated from the highly-doped 
contacts by energy barriers in the access regions due to 
a poor n-type doping caused by masking silicon nitride 
spacers next to the gate (Fig.Q})) . Here we focus on short 
devices with 60 nm gate length where the subthreshold 
regime exhibits specific conductance peaks due to iso- 
lated dopants in the channel. All the measurements pre- 
sented here were gathered on the same sample, but sim- 
ilar features have been observed in other samples of the 
same length (the width, 385 nm here, is not important be- 



cause of a corner-effect [Hi, LL2J). The experiments were 



carried out in a dilution fridge at a base temperature of 
100 mK. The differential conductance was measured with 
a lock-in technique using a 50 /iV amplitude ac voltage. 

The energy barrier between source and drain, mea- 
sured by fitting the thermally activated current at tem- 
peratures above 150 K flol |. is shown on Fig.QJ:. At low 
gate voltage, when the middle of the channel is the high- 
est point of the potential profile, the barrier height de- 
creases linearly with a coupling a — dEb / edV g = 0.68. 
An extrapolation of this linear part to zero barrier shows 
that the conduction band edge in the channel reaches the 
Fermi level at 400 mV gate voltage, value above which the 
low temperature conductance is expected to rise from 
zero. However, several resonances (al, bl, a2, b2) are 
visible at energies below the conduction band edge (inset 
of Fig. [Ik). They come from isolated dopants in the chan- 
nel and will be analyzed in details throughout the paper. 
Fig. [l]; also show that residual barriers remain at larger 
gate voltage. Located in the access regions, they create 
a potential well in the conduction band of the channel, 
with Coulomb blockade at low temperature as revealed 
by the conductance peaks cl and c2. 

When magnetic field is applied (parallel to the current 
here), each peak shifts linearly with field (Fig. [2k ) and 
the slopes are consistent with the Zeeman shift of spins 
±1/2 having a g- factor of 2 as expected for electrons in 
silicon (using the conversion factor a of each peak as 
discussed later). Peaks labelled with a 1 (2) correspond 
to tunneling of spin up (down) electrons under magnetic 
field. Since the first visible resonances are always of spin 
up type (also observed in other samples), we conclude 



that the first peaks (al and bl) correspond to tunneling 
of a first electron (N = — > 1) on two different impuri- 
ties (a and b), and the second set of peaks (a2 and b2) 
correspond to the second charge state (N = 1 — > 2) of 
each impurity (see Fig. [5b). Note that bl can not be the 
second electron of al because a two-electron ground state 
must be a singlet with opposite spins Peaks cl and 
c2 above the band edge are also respectively spin up and 
down under magnetic field. They correspond to the first 
and second electrons of a larger dot (called c) defined 
electrostatically by the gate in the conduction band pro- 
file (Fig.[2j:). In this sample, the higher charge states of 
dot c (N > 2) are not visible probably due to the pro- 
gressive loss of Coulomb blockade above 460 mV where 
the confining barriers disappear. In contrast, in all sam- 
ples, the impurities like a and b always have exactly two 
charge states (N < 2) as expected for dopants because 
their single positive charge can bind no more than two 
electrons and with a very small second binding energy 
[3]. In addition, the measured Zeeman shifts are exper- 
imentally insensitive to the field direction as expected 
for the 3D Coulomb potential of dopants as opposed to 
quantum dots in two-dimensional electron gases (2DEG) . 

Isolated n-type dopants in the p-type channel may 
come from residual impurities in the SOI layer or from 
ion implantation of the contacts. Those contributing to 
the conductance should be close to the gate because of 
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FIG. 2: (color online), (a) Differential conductance plotted in 
color scale (optimized for each peak) showing the peak shifts 
under magnetic field and the g — 2, S — 1/2 Zeeman shifts 
(solid lines) for spin up or down (charging from to 1 or 1 to 2 
electrons), (b) Spin polarization of successive resonances and 
interpretation in term of two dopants a and b, and one larger 
dot c. (c) Conduction band profile with attracting Coulomb 
potential of isolated dopants, for an applied gate voltage cor- 
responding to resonant tunneling through a single dopant in 
the channel. 
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the strong band-bending and along the edges of the wire 
because this band-bending is stronger due to the corner- 
effect (111 ], This strong restriction leads to a silicon vol- 
ume containing less than one residual donor (considering 
that p-type silicon should contain less donors than accep- 
tors). The observed isolated donors must therefore arise 
from the implantation process. 

These isolated dopants are atomic-like systems with 
one or two electrons bound by Coulomb attraction to a 
central ion and can be characterized by their first and sec- 
ond binding energies. These values are known for dopants 
in bulk silicon, but we show in the following that we can 
also measure them in gated nanostructures. Binding en- 
ergies are related to the gate voltage spacing between the 
conductance peaks and the conduction band edge located 
at 400 mV (Fig.QJ;). However the energy change on the 
dopant site is a fraction a — C g / {C g +C s +C d ) of the gate 
voltage spacing, with C s , Cd, and C g being the capaci- 
tances between the dopant and respectively the source, 
drain, and gate electrodes. This ratio is determined from 
the conductance diagram versus gate and bias voltages 
(Fig.[3jD) where positive and negative slopes of conduct- 
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FIG. 3: (color online), (a) Sub-threshold conductance show- 
ing resonant tunneling through dopants a and b. They have a 
similar peak spacing between first and second electron states. 
Inset: peak bl increasing for decreasing temperatures, show- 
ing the quantum regime of Coulomb blockade, (b) Stability 
diagram, i.e. differential conductance plotted in color scale 
versus gate and bias voltages. Excited states (ES) of dopant 
b are indicated. Other weak parallel lines inside the conduct- 
ing regions result from local density of state fluctuations in 
the contacts. 



ing sectors equal C g j (C g + Cd) and C g /C s (the drain is at 
ground). We get C s = 0.15 C g and C d = 0.15 C g for the 
first peaks (al or bl), and C s = 0.27 C g and Cd = 0.31 C g 
for the second peaks (a2 or b2) probably due to thinner 
barriers at higher gate voltage. The dopants are there- 
fore well centered in the channel (C s ~ Cd) and mainly 
coupled to the gate (C s <C C g ). The resulting values 
a = 0.77 (al or bl) and a — 0.63 (a2 or b2) are consis- 
tent with the coupling factor a — 0.68 determined pre- 
viously from the gate voltage dependence of the barrier 
height in the channel below 300 mV. This result validates 
our assumption that the conduction band in the middle 
of the channel evolves with the same coupling below and 
above 300 mV. This point is important since we use it 
to conclude that the band edge crosses the Fermi level 
at 400 mV (Fig.[TJ;) with two dopants below (peaks al to 
b2) and a larger dot above (peaks cl and c2). In the 
following we use a = 0.7 as an averaged value for the 
conversion factor. 

From the 84 mV (74 mV) spacing between peak al (bl) 
and the band edge, we get 59meV (52meV) for the first 
binding energy of dopant a (b). Since this D° charge 
state is neutral, its binding energy is not expected to 
be modified by the capacitive coupling with the nearby 
electrodes. For technological reasons, both arsenic and 
phosphorus dopants are used to make the contacts. Our 
measured values correspond better to the 54meV bulk 
value of arsenic than to the 46meV of phosphorus [HI]. 
The small difference between a and b could be related 
to fluctuations in the band edge energy in the channel, 
induced for example by disorder at the oxide interface. 

The binding e nerg y of the second electron (2 meV only 
for bulk dopants [l6|) is usually very small in this atomic- 
like system because of a strong repulsive Coulomb inter- 
action called charging energy (52 meV for bulk arsenic 
dopants). In our case however, since the D~ state is 
charged, the capacitive coupling with the nearby elec- 
trodes is expected to reduce the charging energy to a 
value e 2 /C where C = C g + C s + C d . The 46 mV (48 mV) 
spacing between peaks al and a2 (bl and b2) indicated on 
Fig. [2k- corresponds indeed to a charging energy between 
the one-electron and the two-electrons bound states of 
32meV (34meV) lower than the 52meV of the bulk. 
The deduced dopant/gate capacitance C g — 3.5 x 10 -18 F 
(3.3 x 10 _18 F) corresponds to a sphere with 2.7nm ra- 
dius in bulk silicon or near a metallic plane covered with 
a thin layer of silicon dioxide. This radius compares very 
well with the 2.5 nm Bohr radius of a neutral D° donor 
in bulk silicon or with the 3 nm Bohr radius estimated for 
the negatively charged D~ state [17|. As a conclusion, 
the second binding energy of the dopant+gate system 
is significantly increased by the capacitive coupling, and 
this effect could provide a longer life-time for the two- 
electron state involved in the read-out scheme of dopant- 
based quantum gates. 

We now turn to the excited state (ES) spectroscopy of 
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FIG. 4: (color online), (a)-(b) Current steps versus gate volt- 
age from to 10 mV bias. The ground state (GS) and an ex- 
cited state (ES) are visible for bl and b2. (c)-(e) Color-plot 
of the current step al at different magnetic fields showing the 
Zeeman splitting of spin up and spin down contributions, (c) 
single step at T, (d) splitting at 5 mV bias between and 
10 T, (e) double step at 10 T with a fully polarized spin up 
current step. 



the two charge states. In bulk silicon the six-fold valley 
degeneracy of the conduction band is lifted on a dopant 
site by the tetragonal symmetry and the first electron 
ground state (GS) is split into one non-degenerate GS 
and a few ES all around 12meV above the GS for phos- 
phorus atoms (22meV for arsenic) [15]. An ES of bl 
is visible at 5meV on Fig. 2k as a second current step 
on top of the GS one. This lower level spacing than 
expected may originate from the strong band bending 
and Stark effect near the gate [l8| which are inherent 
to these nano-structures. This interesting point will be 
the object of further investigations. Regarding the two- 
electrons state, no ES would be observed for a bulk D~ 
ion since any excitation lies in the continuum (l4^ . Here 
however this charged state is more strongly bound due 
to the capacitive coupling with the nearby electrodes and 
several ES can exist. An ES of b2 is indeed clearly vis- 
ible at 3meV (Fig.2t>) and has a larger amplitude than 
the singlet GS in agreement with the triplet degeneracy 
expected for this first ES with two electrons. 

These level spacings of several meV are quite large as 
compared to 2DEG quantum dots and explain why the 
quantum regime of Coulomb blockade [l9| is observed up 
to 15 K with conductance peaks increasing for lowering 



temperature (inset of Fig. [3k, note that at very low tem- 
perature the peak width is limited by the escape rate 
to source and drain) . Another consequence of the large 
level spacing is the absence of level crossing under mag- 
netic field on Fig. [2k because the 1.2 meV Zeeman split- 
ting at 10 T is still smaller than the level spacing. This 
point is also illustrated on Fig.|4]c, d, e showing the cur- 
rent through the state al for a field sweep from to 10 T. 
The single current step at zero field splits into two steps, 
the left one corresponding to a flow of spin up electrons 
and the right one to the sum of spin up and down. The 
Zeeman splitting at 10 T is given by the bias voltage at 
which the two spin states contribute and we get 1.2 mV 
as expected for a spin 1/2 electron in silicon, but here 
the electron is also bound on a dopant. The dotted lines 
on Fig. mark the region of the Coulomb-blocked spin 
down current. The current in the left region is fully po- 
larized with spin up electrons and the dopant acts there 
as a spin filter. A similar magnetic field behavior has 
been observed in GaAs/AlGaAs quantum dots [20] but 
with a much lower g-factor than in silicon and with a 
possible contribution of orbital effects absent in case of 
dopants. 

In conclusion, we used silicon nanowires to investigate 
the spectrum of single dopants in gated nanostructures. 
These dopants form attracting Coulomb potentials and 
induce resonances in the sub-threshold conductance due 
to electron tunneling through their discrete energy lev- 
els. We observed the D° and D~ charge states, measured 
their binding energy and the capacitive coupling with 
nearby electrodes in case of charged states. By means 
of transport spectroscopy versus bias voltage and mag- 
netic field, we measured the quantum level spacing and 
Zeeman spin splitting. 

We thank F. Koppens and M. Sanquer for useful dis- 
cussions. 
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